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ABSTRACT 



We present a study of non-leptonic B-meson decays, partial widths and CP 
asymmetries, and discuss the possibihties of determining the penguin contri- 
butions through less model-dependent phenomenological analyses. In the last 
section we employ a specific model to give more definite numerical predictions. 



1 Introduction and summary of content 



In non-leptonic decays of B-mesons the penguin contributions of the effective Ai? = 1 
hamiltonian become non-neghgible when the tree contributions are suppressed and they 
are even dominant in certain cases. Moreover for those B decays which could give rise to 
observable CP asymmetries, the penguin contributions might cause large deviations from 
the tree predictions even when numerically much smaller than the latter. In such 
cases a correct estimate of the penguin effects is therefore crucial. We present in this note 
a study of such decays. We will discuss the role of the penguin amplitudes both in the 
case where they are not negligible in the computation of the widths and in the case where, 
although not important to evaluate the decay widths, they are relevant to the analysis of 
CP asymmetries. At least in principle, the experimental informations on the first set of 
decays (essentially those corresponding to the Cabibbo suppressed b— > uus transitions) 
might provide bounds, constraints or even a direct access to those parameters defining 
the penguin contributions to the CP violating asymmetries. We shall stress as much as 
possible procedures which are less model-dependent. 

The calculations make use of the recently developed two-loop effective Ai? = 1 hamil- 
tonian 1^ and of the factorization approximation. The form of the amplitudes are collected 
in Tables I, II. They depend on the Cabibbo-Kobayashi-Maskawa matrix elements ^j, on 
combinations Oj of the Wilson coefficients of the effective hamiltonian and on the hadronic 
matrix elements of the weak currents. In view of the inevitable theoretical uncertainties 
we suggest the experimentalists to mostly rely directly on the parametrized expressions 
in Tables I and II to fit the theoretical parameters from the ratios of the partial widths 
which are or will eventually become available (we give some example). 

The asymmetries for decays of neutral B-mesons into states with CP self-conjugate 
particle composition are given in Tables III. The usual asymmetry parameter (neglecting 
strong phases) Im{qA/pA) depends on the mixing phase ipM and on the amplitudes A-re"^^ 
and Ape^'^'' for tree and penguin respectively. The angles ipx, <^p, and the expressions 
for the ratios Ap/At are reported in Tables III. We recall that in some cases the rate 
asymmetry is reduced with respect to the fundamental asymmetry. 

As mentioned above. Tables III, in conjunction with the previous tables I and II, can 
be used to extract, in a way which is as much as possible model independent, the relevant 
parameters, by forming suitable ratios of partial rates and estimating or bounding the 
asymmetries from experimental data directly. We illustrate such a procedure with three 
examples, but various possibilities may appear and the procedures finally chosen will 
presumably depend on the kind of experimental data that will become available at the 
time of the analysis. 

The full predictivity of the model that had been developed in ref. will be exploited 
in the last part of this note to get at the numbers presented in Tables IV. We are mainly 
interested in obtaining an order-of-magnitude estimate of the rates for the processes con- 
sidered in section 2, just to have a rough idea of the rareness of the considered decay 
modes and of their observability at future dedicated facilities. Indeed, as explained above 
and shown in detail in section 2, we think that, maybe, the best use of the these channels 
could consist, via the study of appropriate ratios of rates, in constraining the parame- 
ters entering the theoretical analysis via the factorization hypothesis. These predictions 
should be taken with they are subject to the details of the particular model and 



1 



to the present inputs. The theoretical errors are difficult to estimate. In some case their 
instability with respect to input parameters (for instance Nc) suggests that they should 
be taken with great caution. 

A very large literature already exists on B-meson decays. We shall refer to a recent 
book edited by Stone Q], and for CP violation to the books edited by Wolfenstein and 
by C. Jarlskog 0. 



2 Parametrization of the amplitudes 



The two-loop effective hamiltonian for AB = 1 transitions can be written as follows 



H. 



NL 



V2 



- Vt,v:,Y.^iO, 



i=3 



+ h.o. 



(2.1) 



where q = d,s, are the Wilson coefficients at the scale fi ~ mb, O^'^, 02'^ are the 
operators: 



O" = iub)v-Aiqu)v-A 
O2 = {qh)v^A{uu)v-A 



01 = {ob)v-A{qo)v-A 

02 = {qb)v-Aico)v~A 



(2.2) 



Oi (i = 3, . . . 6) are the so-called penguin operators 

O3 = {qb)v-AY^{qq')v-A 
Oi = ^{q'h)v-A{qq')v-A 

q' 

O5 = {qb)v-Aj2i^'(l')v+A 

<?' 

Oe = -2 5:(g'(l-75)6)(g(l + 75)g') 



(2.3) 



In the previous formulae q = d,s and {ub)v-A = u-y^{l — 75)^. The values of the Wilson 
coefficients, for nib = 4:.8GeV, rritop = 150 GeV and = 250 MeV are P 



ci = 1.133, C2 = -0.291, C3 = 0.015 
C4 = -0.034, C5 = 0.010, C6 = -0.042 



(2.4) 



The hamiltonian ^Tl] allows to study non leptonic decays of B mesons. In a previous 
paper [^, using the factorization approach, we focused on processes where only the op- 
erators Oi and O2 were relevant, at least for what concerns the rates. Notice that the 
penguin operators have Wilson coefficients very small as compared to those of Oi and 
O2, but they become important when the tree diagram is strongly Cabibbo suppressed or 
when the decay can only occur through penguin diagrams. Here we will be interested in 
such processes. Moreover we will study the influence of the penguins on CP asymmetries 
for neutral B mesons. 
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Let us examine first the Cabibbo suppressed decays; in calculating the rates it turns 
out that the penguin contribution is relevant in the flavour processes b — > uus: in this 
case the tree CKM factor is VubV*^, of the order (A is the Cabibbo angle), while the 
penguin has the factor VtbV*g, of the order A^, which compensates the smallness of the 
Wilson coefficients. In Table 1 we give the amplitudes for all these processes, calculated 
using factorization and the following parametrization of the hadronic matrix elements 

< P{p')\V^\B{p) >= [{p + p'r + ^ , ^ glFi(g^) - ^ , "" q'Foiq') (2.5) 



<y(e,p')l ( - A^)\B{p) >- 
+ i{MB + Mv) 



Mb + Mv 



t -q 



— I 



+ le ■ q- 



{Mb + My 

2Mv 



(P + P')/.- 



A2{q') 



where P is a pseudoscalar meson and V a vector meson, q = p — p' 

Mv 



and Fi(0) = Fo(0). We take also: 



2Mi 



Mpt Mv + Mr 



V 



2Mi 



V 



(2.6) 



(2.7) 



< 0|A^|P(p) > = -ifpp" 
<0\V''\V{p,e)> = fvMve^ 



The coefficients a,- 



i = 1, . . .6) are given by 

^21-1 



, C2i C2i-l . 

C2i-i + = C2i + -j^ z = 1, 2, 3 



(2.8) 



(2.9) 



where in the naive factorization Nc is the number of colors, but due to the uncertainties 
of the approach it is treated as a parameter to be determined by the data. For example 
in D decays the value Nf, = oo agrees with the experiments. In Table II we give the 
amplitudes in the factorization approximation for the pure penguin decays of Bg and 
B~ mesons. 

A possible experimental determination of the coefficients could be done through 
ratios of widths, such that most of the hadronic uncertainties cancel out. As an example 
one could measure 04 from the ratio: 



R 



r(po 



TT 



K* 



r(PO ^ TT+p- 



IKfeV;>i| 



(2.10) 



We remark that the width for B^ — > 7r^K*~ depends on the relative phase between tree 
and penguin contributions; in general one should expect a strong and a weak phase, that 
in this case is 7. As for the strong phases, they are expected to decrease as the inverse 
heavy meson mass and therefore we neglect them. 
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3 CP asymmetries 



The penguin amplitude, through interference with the tree amphtude, contributes to CP 
asymmetries in B decays . Neglecting possible strong phases, we write the amplitude of 
the process as a sum of a penguin and a tree contribution 

A = Are'^P^ + Ape'"^"- (3.1) 

where v^t,p are the weak phases. The time dependent asymmetry of a neutral B meson 
{B^ or Bg) decaying into a CP-eigenstate can be written as 

a{t) = -Im sin(Amt) (3.2) 

where 

A=</|5°> A=</|5°> (3.3) 

q/pis the mixing and Am the mass difference of the two neutral B mesons. If Ap/Ax << 1 
we can make an expansion, obtaining: 

/m(-4] = -sm2{(pM + ^t) -'2^cos2{(pM + ^T)sm{(pp - (fr) (3.4) 
\pAJ At 

where is the mixing phase, defined as 

- = e-2*^^ (3.5) 
P 

If ifT = (fp, the penguin does not influence the asymmetry, as it happens in B^ ipK^: 
but in general one has a phase difference. Notice that the relative contribution of the 
penguin in the asymmetry depends on the values of the angles v^a/, V^t and ipp: even 
when Ap/At << 1 it could be quite sizeable (or dominating). In Table III we give, 
for many B^ and Bg decays into states with CP self-conjugate particle composition, the 
values of the phases (pT and ipp\ we give also the branching ratios as in ref.[[7|, and the 
ratio Ap/At as a function of the coefficients a^, assuming factorization. The angle a, /3, 
7 are the usual ones, and (3' is defined as: 



[5' = arg 



(3.6) 



This angle is very small, | sin2/3'| < 0.06. The knowledge of the values of the coefficients 
allows to evaluate the ratio Ap/At'-, in ref. this one is extracted, for the decay 
5° 7r+7r~, from the ratio of the widths of t:^K~ to 7r"'"7r~. As discussed in the previous 
section the advantage of making ratios is the cancellation of some hadronic uncertainties. 
As examples we will consider the decays Bs —>■ D^Dj, Bg p^Kg and Bg — > K^K~ . 
Other examples can be constructed using our tables. The final choices will depend on the 
relative quality of the experimental information that will become available. 

The amplitude for the process Bg DfD~ is 

A = Atq'^' + Ap (3.7) 
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where 
At 

Ap 

and qq'_ 



G 

-^\VcbV*,\ai < Dt\ch4Bs >< Dj\sc^\0 > 
V2 



Gp ( 



M2 



'(mfc - mc){mc + 
(77^(1 — 75)5'. From we obtain: 

_ sin(2/5') + 2xsin/3' 2/3' 



Im 



1 + + 2x cos /5' 1 + X 



(3.9) 



where we have used the smallness of (3' {(3' < 0.03). The asymmetry is given by |3.2| . The 
value of a; = Ap/Ap can be extracted from the ratio of the widths 



R 



V{B, ^ Dtp-) ^ ([dXWd^ 

T{B,^DtTi-) [fj \p^\ 
2 



(M 



D, 



D„ 



fn J IPttI 



\MB^+MnJ'-My 



i{v ■ v') 



H 2 



(3.10) 



(Ms, +Mz5j2-M2 J [i{vv")_ 

where ^ is the Isgur-Wise form factor, v ■ v' = MpJ {'iMp,^) and v ■ v" = (M^^ + M|,^ — 
M^) / {2Mb^MdJ . We work in the infinity quark mass hmit. 

An interesting decay for the measure of 7 is 5^ ^ P^K^, whose amphtude can be 
written as 

A = Ape-'^ + Ape''^ 



(3.11) 



At = 



A, 



\VubV*a\a2 < Ks\db^\Bs >< p°\uu-\0 > 



Gp 
G 

^\Vtiy*^\a4 < Ks\dh_\Bs >< p°|mu_|0 > 



V2' 



For the CP violating asymmetry for the decay B^ —>■ p^Kg one has 



^pA 

For Bg —>■ p~K^ we have 



^ (qA\ -sin(27) + 2xsin(/5-7) +x2sin(2/?) 
Im — - = 



A 



1 + x2 + 2xcos(/? + 7) 



= v/2('^e-^^-Ape^/' 



V 02 



The ratio of the two widths gives 



[l + a;^ + 2x cos(/9 + 7)] 



(3.12) 



(3.13) 



(3.14) 



(3.15) 



T{Bg p-K+) 2 [1 + alx^/al - Ixa^ja^ cos{/3 + 7)] 

Therefore for each given value of (3 and R we can extract x = Ap/Ap as a function of 7; 
from |3.13| we can after that obtain 7, once the asymmetry, given by p.2| , is measured. The 
ratio a^jax^ with its sign, is obtained by experimental results as for example BRiB' — > 
D^T^-)IBR{B^ D+TT-); for instance in we find 02/01 ~ +0.27. 

Another decay potentially interesting to determine 7 is 1?^ K^K" . The amplitude 

is: 

A = Ape-'^ + Ap (3.16) 



At = ^\VubV*^\ai< K+\ub^\Bs>< K-\su^\0> 
v2 



Ap = %\Vt,V, 



V2' 



tbVtsl a4 



< K+\ub_\Bs >< K-\su_\0 :f3.17) 



Notice that in this case Ap and At are of the same order. The asymmetry is proportional 
to: 



Im 



'qa: 



sin(27) + 2xsin7 



1 + x"^ + 2x cos 7 

Similarly to the previous case one can extract the value of x from the following ratio 



(3.18) 



R 



T{Bs ^ K+K- 



where 



'JkV A^[l + a:^ + 2xcos7] 
J^J [1 + x^y^X^ - 2xyX^ cos(/3 + 7)] 

tan (3 



(3.19) 



(3.20) 



tan (3 cos 7 + sin 7 

As for the ratio |3.15| knowing /3, R and the asymmetry p.l8| one can extract the value of 
7- 



4 Numerical results 

So far we have not exploited the model developed in ref. to attempt at more complete 
(and also more risky) numerical estimates. We will do this in the present section and 
in Tables IV. The numerical evaluation of the rates requires the knowledge of the form 
factors. Moreover one has to chose a value for the parameter A^^ initially introduced as 
giving the number of colors, but usually assumed as a phenomenological parameter to 
account for the uncertainties of the model. The analysis done in seems to indicate 
A^c ~ 2 as an effective value, but the problem is open and the theoretical approach still 
to be clarified. We notice that especially 02, 03 and 05 are strongly dependent on Nc and 
the rates dominated by these coefficients can have large variations. 

The results obtained with the model [§] and Nc = 2 are reported in Table IVa and IVb, 
respectively for the decay of Table I and II. In Table IVa we have simply marked with a star 
(*) the strongly Nc dependent decays, while in Table IVb we have also quoted the values 
with Ac = 3 and A^c = 00, because of the large variations of the branching ratio with Ac. 
We notice that for all the decays of Table I the widths have an interference term between 
penguin and tree depending on cos 7, which at present is unconstrained. Therefore in 
Table IVa we present the branching ratios for the two extreme values cos 7 = ±1. 

Experimentally only upper limits exist at present for the processes in Table IV, and our 
values are below these data. The most stringent bounds as compared with our results are 
H for Br{B- (pK' < 8 x 10"^) and Br{B^ tt+K' < 9 x 10"^). The measure of the 
rate B~ 4>K" , a pure penguin process, will allow to extract the value of [03 + 04 + 051, 
which is quite sensitive to A^c (see Table IVb); the rate for B^ n~^K~ depends on the 
value of cos 7 (see Table IVa), therefore a measure of this width will allow a determination 
of 7, once known the ratio of penguin to tree contribution, which can for example be 
measured in the way suggested in [Q, from the ratio B^ — >■ tt+tt" to B^ —>■ tt^K~ . The 
measure of the rate for B^ —>■ tt^K~ could be an alternative way to determine 7 with 
respect to the asymmetry in Bg decays, which presents an experimental challenge and 
needs time-dependent measurements. 
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Tables Captions 



Table I Amplitudes for non-leptonic decays with tree level contribution strongly Cabibbo 
suppressed (KfeKTs)- 

Table II Amplitudes for non-leptonic decays that take place only through penguin con- 
tributions. 

Table III and 5° decays into states with CP self-conjugate particle pairs. The tree 
level contribution dominates the width. Tree and penguin weak phases ifT and ipp 
are given together with the rate of the respective amplitudes. 

Table IVa Branching ratios for non-leptonic decays with tree level contribution strongly 
Cabibbo suppressed {VubV*^) are given for two values of the weak angle 7. The 
contribution is the same in the two cases when there is no penguin contribution to the 
amplitude. We take iKbl = 0.03, = \Vts\ = 0.042 and tb = tb, = 14 x IQ-^^s, 
rrtfe = 5000 MeV, rris = 150 MeV, rrid = 10 MeV and = 5 MeV. Furthermore 
we take = We fixed the value = 2, signaling with (*) the cases where 
there is a strong Nc dependence of the result. 

TablelVb Branching ratios within the model of |^ for non-leptonic decays that take 
place only through penguin contributions. As an example of the variation of the 
result with Nc we have listed the three cases Nc = 2, 3, 00. We take Vtd = 0.01. 
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